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The ability to control spins in semiconductors is important in a variety of fields including 
spintronics and quantum information processing. Due to the potentially fast dephasing times 
of spins in the solid state [U [21 E], spin control operating on the picosecond or faster timescale 
may be necessary Such speeds, which are not possible to attain with standard electron spin 
resonance (ESR) techniques based on microwave sources, can be attained with broadband 
optical pulses. One promising ultrafast technique utilizes single broadband pulses detuned 
from resonance in a three-level A system [3]. This attractive technique is robust against 
optical pulse imperfections and does not require a fixed optical reference phase. Here we 
demonstrate the principle of coherent manipulation of spins theoretically and experimentally. 
Using this technique, donor-bound electron spin rotations with single-pulse areas exceeding 
7r/4 and two-pulses areas exceeding tt/2 are demonstrated. We believe the maximum pulse 
areas attained do not reflect a fundamental limit of the technique and larger pulse areas could 
be achieved in other material systems. This technique has applications from basic solid-state 
ESR spectroscopy to arbitrary single-qubit rotations [U E] and bang-bang control [B] for 
quantum computation. 

Ultrafast optical techniques have previously been used to study semiconductor spins. 
Faraday rotation and differential transmission can be used to passively study spin dynamics 
on the picosecond and femtosecond timescale JTJ [8]. Active techniques to date have utilized 
either the optical Stark effect [9] or Raman transitions through the resonant excitation of 
an optically active state [21 QUI [III H2]- In contrast to previously demonstrated Raman 
techniques, in this paper we present a new stimulated Raman transition technique which is 
based on a single pulse far-detuned from the optical transition. By working off-resonance, 
decoherence due to real population in the excited state is eliminated without the strict 
requirement of 2ir area pulses [5] . We first introduce the theoretical basis for this technique 
before presenting the experimental results. 

We consider the general A-type system with multiple excited states as depicted in Fig. [TJ 
The lower states are denoted |1) and |2) and in the neutral donor system consist of the 
spin-up and spin-down states of the bound electron. These states are coupled via optical 
dipole transitions to the (n — 2) neutral donor-bound exciton states labeled \k). To see how 
the spin can be rotated via a single optical pulse, consider the n-level Hamiltonian in the 
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in which u;l is the Zeeman splitting of the lower states and A*, is the detuning of the applied 
pulse from the |1) «-> \k) transition. The Rabi frequency Qki{^k2) is the product of the 
dip ole matrix element for the |1) |fe)(|2) <-> |&)) transition and the time-dependent electric 
field amplitude E 1 . 

This many-level system can be approximated as a two-level spin system by the adiabatic 
elimination of the upper states, which is valid when the detunings A& are much larger than 
other rates in the system [13] . The effective two-level Hamiltonian is given by 
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Population is thus coherently transferred from one lower state to the other at the effective 
Rabi frequency Q e g(t). Note that the optical phase of the pulse is no longer present in 
this two-level reduction, eliminating the requirement of a fixed optical reference phase as 
a clock signal. If Q c s{t) ^ and \Qi\ = jf^l) the rotation axis will be perpendicular to 
the magnetic field, and full ir rotations with a single pulse are possible. This condition may 
often be met by controlling the polarization of the pulse due to the selection rules for the 
|1) «-> I A;) and |2) <-> \k) transitions. In material systems in which perpendicular rotations are 
not possible, large area rotations can still be achieved using multiple pulses by controlling the 
pulse arrival times over multiple Larmor periods. For a single-pulse rotation, the phase of the 
rotation is determined by the phase difference between frequency components separated by 
the Zeeman frequency within the pulse spectrum and thus is determined by the pulse arrival 



time [I]. We present the simplified two- level approximation as an intuitive description of the 
Raman-rotation technique. However, as will be seen below, in a more realistic three-level 
density matrix model which includes excited state relaxation, high fidelity rotations can still 
be obtained in a non-adiabatic regime. 

One such A-system with multiple excited states that is found in all semiconductors is 
the neutral donor-bound exciton system. For the experimental demonstration of the Ra- 
man technique, we focused on an ensemble measurement of electrons bound to donors in 
bulk GaAs. At liquid helium temperatures the donor electron is bound to the donor impu- 
rity creating a neutral donor (D°). The D° complex is an attractive potential for excitons 
(electron-hole pairs) and the resulting neutral donor-bound exciton (D°X) consists of the 
impurity atom, two bound electrons in a spin-singlet state, and a bound hole. The two 
D° spin states and multiple D°X states are connected by strong, optical transitions [H] and 
form the lower and excited states of our n-level A-type system as shown in Fig. [TJ 

In the first experiment we demonstrate population transfer between states |1) and |2) with 
a single pulse in a 7 T magnetic field. The experiment consisted of three steps: initialization 
of the spin population into state |1), fast-pulse spin transfer, and population readout of state 
1 2). To initialize the spin state, a continuous- wave field was applied on resonance with the 
|2) «-> |3) transition for 10 /xs. This transition is the brightest transition in the D°X spectrum 
shown in Fig. [2^,. At the end of the optical pumping pulse, the state |1) population was 0.94. 
After a 2 /xs delay, which was short compared to the longitudinal relaxation between the 
lower two states [T5], a 2 ps pulse was applied. A typical pulse sequence is shown in Fig. j2}o. 
The pulse was detuned 1 THz below the lowest D°X transition. After a second 2 /is delay 
the optical pumping pulse was again applied and the population in state |2) was measured 
by monitoring the photoluminescence (PL) emitted from the |3) — > |1) transition at the 
beginning of the pulse. The photoluminescence trace as a function of time for a typical case 
is given in Fig. [2^. The conversion from PL intensity to state |2) population was made by 
measuring the PL intensity after the system was allowed to return to thermal equilibrium. 

As shown above, the relative magnitude of Qi and Q2 in Eq. [3j which can be controlled 
by the pulse polarization, determines the rotation axis and must be equal for rotations 
about an axis perpendicular to the magnetic field. In order to obtain the most efficient 
population transfer possible, the experiment was first performed at a constant fast-pulse 
energy for varying pulse polarizations. The polarization dependence is given in the Fig. [2ji 
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inset. Subsequent measurements were performed at the peak of this polarization curve. In 
a system in which Rabi oscillations are observed, the amplitude of these oscillations could 
be used to determine the final rotation axis angle. 

In the single-pulse experiment, population transfer was measured as a function of pulse 
energy by varying the average power of the pulse train. The results are shown in Fig. |2]d. 
At low pulse energies there is a non-linear increase in population with energy which is 
characteristic of coherent Raman population transfer. At higher energies, however, popu- 
lation transfer saturates at a value of 0.5. This saturation was not expected based on a 
three-level simulation which includes the experimentally reported relaxation rates for the 
D°-D°X system (see the Methods section). In this simulation, a numerical solution of the 
master equation, 

p = -i[H, p) + C{p) (5) 

is obtained, where p is the three-level density matrix, H is the Hamiltonian given in Eq. |6j 
and C(p) is the relaxation super-operator given in Eq. [7j In the full density matrix formalism 
we find that although the evolution is not adiabatic and virtual excitation of the excited state 
occurs, high fidelity rotations are still possible. Applying a 2 ps full-width-half-maximum 
hyperbolic secant pulse with a detuning A 3 =l THz and using the definition of fidelity 
F = (ip\p\ifj), where is the desired quantum state, we found that tt rotations should have 
been possible with a 0.97 fidelity and ir/2 rotations should have been possible with > 0.99 
fidelity in the case that = However, a good fit to the experimental data can be 

obtained using a three-level model which includes a dephasing rate of the excited state (73) 
that is linearly dependent on the pulse energy. Although Rabi oscillations were not observed 
in this first experiment, the non-linear increase in population transfer suggests that small 
coherent rotations are possible at low powers. 

In order to measure the coherence of the small angle rotations, a second experiment with 
two fast pulses was performed. In the double-pulse experiment the single pulse was split into 
two pulses with a variable delay t d (Fig. |3^). As the delay between the two pulses increases, 
the final population in state |2) oscillates with a period equal to the Larmor period r^. 
In Fig. [3}d we plot the population in state \2) after two pulses of energy 10 pJ/cm 2 were 
applied, as a function of td- An oscillation at the Larmor frequency is clearly observed 
verifying coherent population transfer as well as rotation axis control. The observed 42 GHz 
oscillation at 7 T corresponds to a D° electron g-factor of g e = —0.42 which is consistent 
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with previous measurements [IB] . In contrast to the ideal case, the population in state 
|2) never reaches the optically pumped value and perfect destructive interference does not 
occur. The finite population left in state |2) indicates that in addition to coherent population 
transfer, incoherent population transfer occurs. The simultaneous fit of the single (Fig. [2jl) 
and double-pulse data (Fig. [SJo) to the energy-dependent dephasing model indicates a single- 
pulse rotation of 0.9 radians and a double-pulse rotation of 1.8 radians with fidelities 0.85 
and 0.78 respectively. The two-pulse experiment was performed at several powers with a 
linear decrease in the visibility observed with increasing power as shown in Fig. ^jp. 

The saturation of the population transfer in the first experiment and the absence of 
perfect destructive interference in the second experiment indicate that some fast dephasing 
mechanism is occurring in the D°-D°X system. This behavior is not expected from the 
known parameters in the D°X system in the low excitation limit [T]. However the single 
impurity three-level model may not be sufficient to describe the experimental many-exciton 
system. We note that in addition to exciting an appreciable virtual D°X population during 
the applied pulse, virtual free excitons are also excited. As shown in Fig. [2^i, the D°X tran- 
sitions lie right on the tail of the free exciton transitions. Unlike the D°X transitions, the 
free-exciton transition is quite broad (THz) and thus additional real excitation due to the 
picosecond pulse is likely. We expect the incoherent free exciton excitation to be linear in 
pulse-energy until the free-exciton transition has saturated. Once free-excitons are excited, 
exciton-exciton interactions and exciton-electron interactions [TTl [T%] could be the source 
of the observed fast dephasing and a many particle model may be necessary to explain the 
experiment result. Dephasing due to multi-exciton interactions should be less of a factor in 
deeper impurity systems [19] or charged III-V quantum dot systems [2J [TOj [11] which have 
much larger exciton binding energies than the GaAs D°X system. In these systems it may 
be possible to obtain single-pulse large area spin-rotations which would be a valuable tool 
for spin-based quantum information processing. 

Even with the modest population transfer possible in our bulk experimental system, 
simulations indicate that 7r-pulses are still possible if several low energy pulses are applied 
to the system in phase. Again, using the three-level model that assumes a level |3) dephasing 
rate (73) with a linear power dependence, we find that a 7r-rotation is possible with 0.80 
fidelity in as few as 8 pulses (200 ps) with energy densities of 5 /iJ/cm 2 each. As visible in the 
Fig. [4] inset, applying small-angle pulses in succession results in a discrete Rabi oscillation 
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curve. Fig. [4] shows the fidelity of a tt rotation as a function of the number of pulses applied 
in phase. One might expect the fidelity to increase to 1 as more and more weaker pulses are 
used to perform the rotation. However the finite T 2 decoherence time of 1 ns limits the total 
number of pulses that can be used. 1 ns is the experimental inhomogeneous dephasing time 
T 2 * in the system [TJ and the decoherence time could be much longer for a single spin system. 
While the multiple-pulse technique is significantly slower than the single-pulse technique, it 
may prove valuable in systems in which high-area rotations are desired yet low powers are 
necessary. Such large area rotations are necessary for single gates in quantum computation 
and can also aid in the suppression of decoherence [6]. However, for general ESR techniques, 
such as spin-echo, small area pulses are all that is necessary for determining the fundamental 
homogeneous decoherence time T 2 [20J in a material. Thus this Raman fast-pulse technique 
as experimentally demonstrated can immediately be applied for these studies. 

I. METHODS 

A. The D°X system 

The sample studied consisted of a 10 /mi GaAs layer with a donor density of 5 x 
10 13 cm -3 on a 4 /iin Alo.3Gao.7As layer grown by molecular-beam epitaxy on a GaAs sub- 
strate. The sample was mounted strain-free in a magnetic cryostat in a liquid helium bath. 
The magnetic field was parallel to the < 110 > crystallographic axis. The magnetic field 
was perpendicular to the excitation and collection paths. The signal was collected from a 
20 yum spot with an estimated 10 5 donors contributing to the signal. 

In the applied magnetic field, the neutral donor electron splits into two levels denoted by 
the magnetic quantum number m e = ± | as shown in Fig. [TJ The bound electron g-factor g e 
is -0.40 to -0.46 depending on the strength of the magnetic field and the crystal orientation 
with respect to the magnetic field [16] . The excited state D°X energy level structure is much 
more complex and although it is not fully understood a detailed study has been performed 
by Karasyuk et al.[16j. The two electrons in the complex form a spin singlet denoted in 
Fig. [I] by \m s = 0). The energy of the D°X state is thus determined by the spin of the 
bound hole nih = ±| as well as the hole's effective mass orbital angular momentum L. 
In our system we have identified the ground and first excited state of the D°X complex as 
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the \L = l,rrih = — \) and the \L = 0,rrih = — |) states respectively. 

A GaAs photoluminescence spectrum with above-band excitation at 815 nm can be seen 
in Fig. [2^ (blue curve). Only the horizontal polarization is collected to resolve more of 
the D°X transitions and D°X linewidths are instrument resolution limited. Also observed 
in the spectrum are the free exciton transitions, the acceptor bound-exciton (A°X) transi- 
tions, and the D°X two-electron satellite transitions (TES). TES transitions occur when the 
D°X relaxes into a D° excited state. 

High resolution photoluminescence excitation (PLE) spectroscopy as well as the location 
of the TES lines indicate that the primary donor in our sample is silicon with a much 
weaker (< 10 x) concentration of sulfur. PLE spectroscopy on the sample shows that the 
inhomogeneous linewidth is < 10 GHz on the broadest lines pQ. This linewidth is much 
narrower than the 1 THz detuning and thus should not affect the fidelity of the pulse 
rotations. 

B. More on the experimental pulse sequence 

A typical fast-pulse sequence detected on a GHz photodiode is shown in Fig. [2]d in 
which the pump/read-out pulse intensity is 200 times the intensity used in the experiment. 
The optical pumping laser was a Coherent Ti: Sapphire 899-29 continuous- wave (CW) laser 
modulated by an AOM and was polarized parallel to the magnetic field. Depending on the 
experiment the pump intensity ranged from 5 to 15 /iW and the spot size was approximately 
120 /im. The fast pulse was provided by a picosecond mode-locked laser with a repetition 
rate of 80 MHz. Single pulses were picked by an EOM every 10 to 15 /is depending on 
the particular experiment. The EOM extinction ratio ranged from 60-100 and an additional 
2 /as extinction envelope was provided by an AOM (extinction ratio of 1000). The fast-pulse 
polarization was 45 degrees from the magnetic field axis. The fast-pulse laser spot size was 
80 microns. 

C. Three-level theoretical model 

The data in Figs. [2]d and [3]d are fit to a three level density matrix model which solves 
the master equation given in Eq. [5j The three-level Hamiltonian in the rotating frame is 
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given by 
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in which Ti2 (r 2 i = T^e^kr) is the longitudinal relaxation rate from |1) — > |2) (|2) — > |1)), 
2r 3 is the total radiative relaxation from |3), 72 is the transverse relaxation rate between 
|1) and 1 2), and 73 (t) is the level |3) dephasing. For electrons bound to neutral donors 
in GaAs, the relevant relaxation parameters are the excited state radiative relaxation rate 
2r3 = (1 ns) _1 [II], the lower state decoherence time T 2 = 1 ns [lj, and a Zeeman splitting of 
ul = 42 GHz in a 7 T field. In the constant level- 1 3) dephasing model 73 (t) = 10 GHz [T1IT6]. 
Three level simulations show that even in this non-adiabatic regime, D°X population is only 
virtually excited during the pulse duration. Pulse fidelities are high as long as all dephasing 
mechanisms in the system are slow compared to the pulse time. In the intensity dependent 
model the peak value of 73 (t) is given in Fig. [2ji. The model fits the data if pulse energies are 
a factor of 0.8 less than the energies in the experiment. This discrepancy could be explained 
by our incomplete knowledge of the D°X dipole matrix elements. The model includes only 
one excited state with a relaxation rate equal to the experimental total relaxation rate for 
the many level D°X state at zero magnetic field. In reality our fast pulse is interacting with 
many levels in a 7 T field. 
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FIG. 1: ra-level energy diagram with an applied time-dependent electric field. The two ground 
levels are split in energy by %ujl. An electric field E with energy Uojq is applied to the system 
detuned by the energy ft A3. Energy levels corresponding to the donor bound exciton system are 
denoted to the right of the diagram and are explained further in the Methods section. 
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FIG. 2: Description of the single-pulse experiment and results, a) Black line: GaAs photolumi- 
nescence spectrum with above-band excitation at 815 nm. Grey line: Fast pulse spectrum tuned 
1 THz from the lowest D°X transition, b) A typical laser-pulse sequence recorded on a GHz pho- 
todiode. c) A typical PL trace collected from the |3) — * |1) transition as a function of time, d) 
Markers: Experimental population in state |2) as a function of pulse energy density. Data were fit 
to an exponential to determine the initial PL intensity and thus the initial population in state \2). 
Error bars are derived from the standard deviation of the fitting coefficients. A saturation at high 
energies is observed. The arrow marks the pulse energy used for the two- pulse experiment shown 
in Fig[3j Dashed line: Three-level simulation assuming a constant level- 1 3) dephasing rate. Solid 
line: Fit of the data to a three-level simulation using the same parameters as before, except for a 
state |3) dephasing (73) with a linear energy dependence. Dotted line: Peak value of 73 used for 
each pulse energy density in the energy dependent dephasing model. Inset: Population measured 
in state |2) as a function of pulse polarization. 
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FIG. 3: Description of the double-pulse experiment and results, a) Schematic diagram of the 
double-pulse experiment pulse sequence, b) Markers: Population in state \2) after the second 
pulse as a function of T£>. Solid black line: The simulation result that uses a level |3) dephasing 
of 1.6 THz consistent with the energy dependent dephasing model. Dashed line: Three-level 
simulation with a 10 GHz level |3) dephasing rate. Dotted line: Residual population remaining in 
state 1 2) after the optical pumping pulse, c) Pulse energy dependence of the two-pulse visibility. 
The solid line is a linear fit. The visibility is defined as r max ~{ m '" in which I max (Imin) is the 
maximum (minimum) of the two-pulse visibility curve after subtracting the intensity observed 
after optical pumping. 
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FIG. 4: Theoretical calculation of the fidelity of it pulse vs. number of pulses applied in phase. 
Inset: Population in state |2) vs. time as 2 ps pulses are applied in phase. 
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